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ABSTRACT 


Several  types  of  surface  interpolation  techniques  exist  which  are 
defined  over  triangular  or  rectangular  surface  patches.  They  all  have 
proven  their  usefulness  in  Computer  Aided  Design. 

This  paper  describes  a  first  approach  to  a  programming  system  that 
combines  triangular  and  rectangular  interpolation  methods  and  applies 
them  to  a  typical  design  problem  in  shipbuilding,  the  representation 
of  a  ship  hull. 

The  system  is  designed  such  that  it  requires  only  those  data  which  are 
available  from  a  common  shiplines  graph  and  produces  the  output  in  the 
most  general  form,  i.e.  as  a  set  of  points  (X,Y,Z(X,Y))  for  the 
single  patches.  This  output  will  then  be  processed  in  a 
post-processor  fashion  by  some  particular  graphical  or  production 
device . 
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INTRODUCTION 


Computer  Aided  Geometric  Design  (CAGD)  is  a  specialized  field  in 
Computer  Aided  Design  (CAD)  and  focuses  on  the  mathematical 
representation  of  arbitrily  shaped  univariate  curves  and  bivariate 
surfaces.  It  provides  several  interpolation  methods  which  already 
have  been  incorporated  in  the  design  process  in  many  companies. 

Interpolation  over  rectangular  and  triangular  surface  patches  are 
two  important  schemes  used  for  surface  representation.  They  have  been 
used  mainly  to  provide  the  surface  data  of  arbitrarily  shaped  objects 
for  the  purpose  of  stress  analyses  and  aerodynamic  computations. 
However,  the  two  schemes  have  not  been  applied  in  conjunction  with  one 
another  yet  to  solve  design  problems.  Therefore  optimal  design 
features  for  general  purpose  applications  have  not  been  achieved  yet, 
because  previous  design  systems  were  strongly  affected  by  the 
advantages  and  disadvantages  of  the  single  methods. 

One  of  the  objectives  of  the  project  to  be  described  in  this 
paper  was  to  combine  the  advantages  of  these  two  methods  and  to 
eliminate  or  at  least  reduce  the  disadvantages  to  a  minimum.  Thus  the 
efficiency  of  the  rectangular  interpolation  schemes  is  exploited  in 
parts  of  the  surface  (e.g.  mid-section  of  a  ships  hull)  where 
rectangular  patches  can  be  used  without  any  problem.  For  the  more 
complex  shaped  parts  of  the  surface  (e.g  bow  and  stern  of  a  ship)  the 
more  elaborate  but  therefore  more  adaptable  triangular  schemes  are 
employed  (see  Fig.  1) . 
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Fig.  1:  Areas  of  rectangular  and  triangular  interpolation 
on  a  ship  hull  surface 
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LINES  AND  OFFSETS 


The  programming  system  for  the  surface  representation  is  not 
intended  to  be  used  for  interactive  ship  hull  design  or  fairing.  It 
is  meant  for  the  more  production  oriented  phase  of  the  entire  design 
process,  at  which  stage  a  fairly  accurate  lines  graph  (e.g.  ship 
lines  graph)  is  expected  to  be  available. 


1.  PROGRAM  DESCRIPTION 

1.1  Input  Data  for  the  Program  SHIPDS 


0  1 

C°  -  and  C-data  which  can  be  obtained  from  the  ship  lines  graph 
are  the  only  input  to  the  programming  system.  These  are  the 
positional  (Xi,  Yi,  Z(Xi,Yi))  and  the  slope  (DZ  (Xi, Yi) /DX, 
DZ  (Xi,  Yi) /DY)  data  at  each  of  the  vertices  of  the  single  patch  i  (see 
Fig.  2)  . 


patch  i 


DZ(X.,Y.)/DY 


D  Z ( X i ,Yi )/DX 
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9*  WL 


4-CW.L 
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Half  -  Breadth  Plan 


Fig.  2:  C°-  and  C1  -  data  obtained  from  the  ship  lines  graph 


The  patches  can  be  defined  easily  on  the  grid  that  is  used  to 
draw  the  ship  lines.  But  it  does  not  need  to  be  the  same  grid  if 
bigger  patches  are  sufficient  for  the  required  interpolation  accuracy. 

The  principal  input  sequence  in  form  of  data  cards  is  shown  in 
Fig.  3.  Fig.  4  illustrates  this  sequence  in  the  form  of  records  in  a 
data  file  for  the  particular  case  of  harbor  tug  Huntington  / 7 / . 
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LINES  AND  OFFSETS 


NDREC  =  number  of  rectangular  patches  in  a  set 

NDTRIA  =  number  of  triangular  patches  in  a  set 

I  PTYPE  =  patch  type,  3:  triangular  patch,  4:  rectangular  patch 

Nl  NT  =  number  of  intervals  on  the  interpolated  patch 

X,  Y,  Z  =  position  at  each  vertex  of  the  patch 

DZ,  DZ  =  gradient  at  each  vertex  of  the  patch 
DX  DY 

XMI  N,  YMI  N,  ZMI  N  =  lower  limits  for  the  X,  Y,  Z  data 

XMAX,  YMAX,  ZMAX  =  upper  limits  for  the  X,  Y,  Z  data 


Fig.  3:  Principle  sequence  of  input  data  for  the  program  SHI  PDS 
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Fig.  4:  Example  of  an  input  sequence  illustrated  on  the 
harbor  tug  Huntington 
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There  is  no  restriction  imposed  on  the  order  of  patch  definition 
nor  on  the  number  of  patches  in  a  set.  The  variable  NINT  controls  the 
grid  size  on  which  the  interpolation  values  are  computed  or,  in  other 
words,  controls  the  density  of  interpolation  points  on  a  surface 
patch.  The  values  zero  for  the  patch  type  IPTYPE,  number  of  patches 
NPATCH  and  number  of  intervals  on  a  patch  NINT,  signal  the  end  of  the 
input  stream  and  terminate  the  execution  of  the  program  SHIPDS. 


1.2  Output  Data  from  the  Program  SHIPDS  and 
Input  Data  for  the  Program  SHIPLO 


All  patches  are  processed  independently  and  sequentially  as  they 
are  defined  in  the  input  stream  (see  Fig.  3  and  Fig.  4) .  Therefore 
the  sequence  of  output  data  from  the  program  SHIPDS  is  of  similar 
form.  Each  interpolated  patch  is  described  accordingly  by  a  set  of 
triples,  the  coordinates  of  the  interpolated  points  X,Y,Z(X,Y). 
Fig.  5  illustrates  the  principal  arrangement  of  the  output  data  in 
form  of  data  cards.  Fig.  6,  however,  demonstrates  this  arrangement 
for  the  particular  case  of  harbor  tug  Huntington  in  the  form  of 
records  in  the  output  file. 

Here  also  the  values  zero  for  the  patch  type  IPTYPE  and  the  number  of 
interpolated  points  on  the  patch  NDRES  indicate  the  end  of  the  output 
file  from  the  program  SHIPDS. 

The  output  data  from  the  program  SHIPDS  are  also  treated  as  input 
data  to  the  program  SHIPLO,  which  represents  any  kind  of  graphical  or 
production  device. 

It  should  be  mentioned  that  the  interpolation  points  are  computed 
on  a  grid  of  lines  with  constant  values  of  the  free  variables.  This 
property  normally  can  be  exploited  by  the  interface  program  for  the 
plotting  device.  In  case  of  rectangular  patches  the  grid  of  the 
interpolation  points  coincides  with  the  X  =  const.  and  Y  =  const, 
lines  (see  Fig.  7)  .  In  the  case  of  triangular  patches,  however,  the 
points  are  computed  on  lines  with  constant  barycentric  coordinates. 
Those  lines  are  parallel  to  the  edges  of  a  triangular  patch  (see 
Fig.  7)  . 


1.3  Structure  of  the  Programming  System  SHIPDS  -  SHIPLO 


The  system  is  called  a  "two  phase"  programming  system,  since  it 
is  separated  into  two  distinct,  and  to  some  extent  independent,  parts. 

SHIPDS  represents  the  "first  phase"  of  the  system  and  performs 
the  interpolation  computations  over  the  rectangular  or  triangular 
patches.  The  main  objectives  for  the  design  of  this  part  of  the 
system  were  portability  and  flexibility.  This  could  be  accomplished 
by  using  standard  FORTRAN  IV  as  the  most  widely  used  programming 
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Fig.  5:  Principle  sequence  of  output  data  from  the  program  SHI  PDS 
and  sequence  of  input  data  for  the  program  SHI  PLO 
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Fig.  7:  Interpolation  grid  of  a  rectangular  and  triangular 
surface  patch 


language  in  engineering  and  by  following  strictly  the  rules  of  the 
standard. 

Exploiting  the  FUNCTION-  and  SUBROUTINE-f eatures  the  program  SHIPDS 
could  be  designed  in  a  very  modular  structure  (see  Fig.  8)  ,  such  that 
the  rectangular  and  triangular  interrpolation  schemes  as  a  whole  could 
be  replaced  by  other  techniques  easily.  This  also  permits  the  use  of 
those  interpolation  programs  in  other  than  just  the 
SHIPDS-environment . 

The  programs  RPINT  for  rectangular  patch  and  TPINT  for  triangular 
patch  interpolation  themselves  are  of  very  modular  structure.  Thus 
different  techniques  for  interpolating  the  edge  functions  of  a  patch 
or  approximating  tangential  and  cross  boundary  derivatives  along  the 
edges  could  be  substituted  by  better  and/or  more  appropriate  ones. 

The  freedom  of  choice  for  the  different  interpolation  methods, 
however,  is  restricted  by  the  type  of  data  they  need.  Since  most 
techniques  do  not  require  more  than  the  C  -  and  C  -  data  as  described 
in  chapter  1.1  at  the  vertices  of  a  patch  and  since  those  data  all  can 
be  obtained  from  the  ship  lines  graph  (see  Fig.  2),  the  system 
deliberately  was  confined  exclusively  to  those  data. 

The  output  data  from  SHIPDS  and  also  input  data  for  SHIPLO  had  to 
be  provided  in  a  most  general  form,  because  the  type  of  graphical  or 
production  device  is  not  known  in  advance  and  will  vary  from 
application  to  application.  The  essential  information  each  kind  of 
graphical  or  production  device  will  need  is  the  coordinates  of  the 
interpolation  points  on  a  patch  together  with  the  number  of  points 
thereon  and  the  code  for  the  type  of  the  patch  (see  Fig.  5  and 
Fig.  6)  . 
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8:  Structure  of  the  programming  system  SHIPDS  -  SHIPLO 


second  phase 


The  "second  phase"  of  the  programming  system,  represented  by  the 
program  SHIPLO,  is  highly  device  dependent.  It  is  an  interface 
program  that  needs  be  modified  or  completely  rewritten  new 
application  of  the  system.  It,  therefore,  was  written  as  a  single 
program  and  very  closely  adjusted  to  the  output  device. 


2 .  THEORY 

2.1  Interpolation  Schemes  over  Rectangular  Patches 


The  parametric  interpolation  method  over  rectangles  was 
introduced  by  S.  A.  Coons  / 1  / .  Therefore  they  are  also  called  "Coons" 
patches . 

This  method  uses  the  boolean  sum  (P1(5p2)F  that  interpolates  to 

the  values  of  the  corners  and  edges  of  a  rectangular  patch.  It  is 
defined  as  follows: 


(P10P2)F  =  P1F  +  P2F  -  P1P2F 


(1) 


with  the  projectors 
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and  the  cubic  Hermite  functions 


hO  (t)  =  (l-t)f  (  2t+l )  hO  (t)  =  (l-t)f-  t 

hi  (t)  =  tA  (-2t+3)  hi  (t)  =  t'Mt-l),  t€[0,l]. 


(2) 


(3) 


(4) 


(5) 


(6) 
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F(u,0) ,  F(u,l) ,  F  (0 ,  v) ,  F  (1 ,  v) 


(7) 


denote  the  edge  functions  of  the  natch  (see  Fio.  9)  and 
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[DF(l,v)/Pv 

denote  the  cross  boundary  derivatives  alono  the  edges. 


Fiq  9:  Bicubically  blended  Coons  patch 


The  edge  functions  are  also  interpolated  by  cubic  Hermite 
functions  (see  (6))  using  the  C  -  and  C  -data  provided  as  input  for 
each  patch  (see  chap.  1.1): 


F(U,1) 

=  jh0(u) 

hi  (u) 

hO  (u) 

hl(uj] 

F  (0,1) 

F  (1,1) 
F10  (0,1) 
F10(l,l)_ 

Fd,v) 
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hO  (v) 
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F  (1,07 
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F01  (1,0) 
F01  (1,1) 
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(10) 


The  cross  boundary  derivatives  are  linearly  blended  together: 


F01 (u,l)  =  (1-u)  F01  (0,1)  +  u  F01 (1,1) 


(11) 


F10  (1  ,v)  =  (1-v)  F10 (1,0)  +  v  F10 (1 f 1 )  (12) 

The  use  of  the  cubically  blended  Coons  patch  generally  suffers  from 
the  approximation  of  the  twist  partition  in  the  matrix  B: 


m 


(13) 


In  order  to  keep  the  problem  simple,  the  twist  partition  T  is 
approximated  by  zero.  But  it  normally  does  not  result  in  a  sufficient 
surface  interpolation.  One  of  several  different  techniques  for 
approximating  the  twist  terms  is  the  so  called  "Gregory's  Square" 
which  is  a  12-Parameter  scheme  (see  / 2 /  for  more  details)  . 

Since  adjacent  patches  have  the  same  positions  and  gradients  at 
the  patch  vertices  as  obtained  from  the  ship  lines  graph  and  since 
their  common  edge  functions  are  interpolated  by  the  same  cubic  Hermite 
functions  and  the  cross  boundary  derivatives  by  the  same  linear 
blends,  the  surface  represented  by  rectangular  patches  must  be  a 
Cl-surface . 


2.2  Interpolation  Schemes  over  Triangular  Patches 


Barnhill,  Birkhoff  and  Gordon  (BBG)  initiated  in  1973  / 9 /  the 

"triangular"  Coons  patch  by  applying  the  boolean  sum  to  the  three 
projectors  of  the  standard  triangle  (see  Fig.  10) : 

(Pi  ©  Pj)F  =  PiF  +  PjF  -  PiPjF  i,  j  =  l ,2,3  (14) 

i#j 

with 

PIF  =  [hO  ( r )  hi  (r)  hO(r)  hi  (r)|  [f  (0  ,ci) 

F  (l-a,q)  (15) 

FJ  0  (0  , q) 

F10(l-q,q)  and 

r  =  P/(l-P)  .  (15) 

P2F  and  P3F  are  defined  analogously.  This  interpolant  interpolates  to 
the  values  at  the  corners  and  edges  of  the  standard  triangle. 

In  order  to  apply  the  BBG-interpolant  to  an  arbitrary  triangle 
affine  transformations  on  the  standard  triangle  need  to  be  performed. 
Positions  and  directions  of  the  tangential  derivatives  are  preserved 
under  affine  transformation  but  not  the  directions  of  the  cross 
boundary  derivatives.  Therefore  they  need  special  treatment  as  , 
described  in  more1  details  in  / 2 /  and  / 4 / . 

The  "Brown  -  Little"  triangle,  however,  is  a  C4-interpolation 

scheme  for  arbitrary  triangles  that  interpolates  to  the  values  at  the 
corners  and  edges  of  the  triangular  patch: 

(BL)  F  =  A*P1  +  B*P2  +  C*P3  (17) 

with  the  weight  functions 

A  =  b22*b32,D  B  =  bl*b3%  n  C  =  bl  *b2*"/D  and 
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bl  *b3 


2  2 

+  b2  *b3 


(Ifl) 


The  bi,  i=l,2,3  are  the  barycentric  coordinates  at  the  point  V 
with  the  property: 


bl  +  b2  +  b3  =  1  (19) 

The  physical  coordinates  X,  Y  can  be  computed  by: 

3 

X  bi*Xi ,  Y  =£jbi*Yi  .  (20) 

4*1 

The  projectors  PiF  are  defined  as  follows: 

PiF  (X, Y)  =  [ (F (pi)  +  bi*2AR) / [Vi  +  1  -  Vi  +  211  *  DF(pi)/Dni  (21) 


where  AR  denotes  the  area  of  the  projection  of  the  triangular  patch 
onto  the  XY-plane  and  DF(pi)/Dni  the  cross  boundary  derivative  along 
edge  i.  The  Vi's  are  the  patch  vertices  and  the  pi's  are  the 
parameters  which  vary  between  the  values  of  the  vertex  coordinates  , 
which  enclose  the  same  edge.  They  are  the  projections  of  an  arbitrary 
chosen  interpolation  point  V  within  the  patch  onto  the  single  edges 
(see  Fig.  11)  . 

In  the  case  of  the  triangular  patches  the  edge  functions  are  also 
interpolated  by  cubic  Hermite  functions  using  the  C°  and  C^-data  from 

the  graph.  Since  adjacent  patches  need  to  have  at  least  one  edge  in 
common,  the  two  corresponding  vertices  need  to  have  the  same 
positional  and  slope  data  at  those  points,  which  are  computed  from  the 
gradients.  Therefore  the  edge  functions  of  adjacent  triangular  but 
also  adjacent  triangular  and  rectangular  patches  will  be  identical  and 
the  surface  constructed  from  the  trianqular  patches  will  be  C*. 


3.  RESULTS 


3.1  Computational  Problems 


In  the  first  stage  of  designing  the  programming  system 
SHIPDS-SHIPLO,  most  of  the  concern  was  dedicated  to  the  program 
organisation  (see  chap.  1.3)  rather  than  to  the  perfection  of  the 
interpolation  programs.  Therefore  those . interpolation  programs,  which 
were  already  available,  have  been  adopted  accordingly  to  the 

requirements  of  the  program  organisation.  Thus  the  bicubically 
blended  Coons  patch  with  the  zero  twist  partition  (see  (1)  and  (5)) 
for  the  rectangular  patches  and  the  Brown-Little  triangle  (see  (16)  ) 
for  triangular  patches  have  been  implemented. 

The  results  are  shown  in  Fig.  12.  They  are  surely  not  satifactory  yet 
for  most  "real  world"  applications.  The  poor  quality  of  the 
interpolation  is  partly  caused  by  the  inaccuracy  of  the  input  data, 
which  have  been  manually  read  off  the  graph  shown  in  Fig.  1  and 
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Fig.  2. 

The  cubic  Hermite  functions  for  interpolating  the  edge  curves  of 
the  surface  patches  have  been  used  in  nonparametric  form  so  far.  This 
also  creates  problems,  since  their  interpolation  accuracy  decreases 
very  fast  with  the  increase  of  the  slopes  of  the  tangents  at  the  patch 
corners . 

*  The  tangential  and  cross  boundary  derivatives  needed  for  both  the 
triangular  and  rectangular  patches  are  computed  from  the  gradients  at 
the  vertices.  In  case  of  the  rectangular  patches,  in  fact,  they  are 
identical  with  the  gradient  components  and  this  part  of  the  ships  hull 
surface  indeed  is  a  C1-  surface  and  interpolates  correctly  over  its 
entire  region. 

The  triangular  interpolation  scheme  also  interpolates  a  C* 
surface  as  expected.  But  the  local,  (  i.e.  in  the  interior  of  the 
patch)  interpolation  accuracy,  especiallv  in  parts  of  the  surface  of 
high  curvature  as  the  stern  of  the  ship,  is  very  low.  Noncollinear 
edges  of  triangular  patches  along  the  surface  boundaries,  on  the  other 
hand,  create  an  inconsistency  for  the  local  interpolation  along  the 
boundaries.  This  problem  has  not  yet  been  solved  theoratically  and 
results  in  inaccurate  interpolation  along  the  surface  boundaries.  It 
will  depend  on  the  type  of  the  surface  as  well  as  on  the-  particular 
application  wether  the  final  interpolation  will  be  sufficient  or  if 
other  correction  steps  have  to  be  taken  (see  Fig.  12  and  Fig.  13). 

Discontinuities  in  a  smooth  surface  as  they  might  be  intended  and 
occur  in  the  case  of  the  harbor  tug  Huntington  (see  Fig.  1  and 
Fig.  2),  can  successfully  be  treated.  But  they  have  to  be  restricted 
to  follow  exactly  the  edges  of  just  the  rectangular  patches. 


3.2  Graphical  Problems 


As  example  for  the  output  device  for  the  "second  phase"  of  the 
system  a  TEKTRONIX  -  storage  tube  has  been  chosen,  that  was  equipped 
with  a  transformation  package  to  create  3-D  images. 

In  order  to  obtain  a  reasonable  impression  of  a  3-D  image  on  a 
2-D  TEKTRONIX  -  screen  a  regular  texture  on  the  surface  patches  is 
needed.  This  can  be  accomplished  easily  for  rectangular  and 
triangular  patches  separately.  But  since  the  interpolation  points  of 
the  triangular  and  rectangular  patches  are  computed  on  different  grids 
(see  chap.  1.2  and  Fig.  7)  it  is  not  easy  to  achieve  a  completely 
homogeneous  texture  on  this  kind  of  graphical  device  and  thus  avoid 
the  considerable  distortion  of  the  principally  correct  results. 
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z 


Y 


Harbor  tug  Huntington  (3-D  front  view) 


Harbor  tug  Huntington  (side  view) 
g.  12:  First  results  from  the  programming  system  SHIPDS  -  SHIPLO 
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4.  NECESSARY  AND  POSSIBLE  IMPROVEMENTS 


In  order  to  obtain  results  that  are  acceptable  for  'real  world" 
applications  of  the  system  some  improvements  have  to  be  achieved. 
Thus  the  rectangular  interpolation  scheme  with  the  'twist  approximation 
equal  to  zero  has  to  be  replaced  by  a  better  me  t  h  o  d . 

The  triangular  interpolant,  the  Brown-Little  triangle,  is  known  for 
its  poor  approximation  of  the  interior  of  the  patch.  It  also  will  be 
substituted  by  an  acknowledged  better  method;  Because  of  the  modular 
structure  of  the  program  SHfPDS  those  changes  can  be  performed  easily. 

The  nonpara  me  trie  interpolation  used  for  the  edges  of  the 
triangular  and  rectangular  patches  will  be  sufficient  for  many 
problems.  It  is  desirable,  however,  to  implement  a  parametric  version 
in  order  to  be  more  general  and  to  avoid  problems  with  extrem  slope 
values  as  they  occur  at  the  bottom  and  the  stern  of  a  ships  hull. 


As  Fig.  13  can  prove  those  "better"  interpolation  s  c  h  e  me  s  are 

avilable  as  part  of  the  surface  editor  SURFED  /  8/  of  the  CAGD  group 

at  the  University  of  Utah.  The  triangular  and  rectangular  methods  can 
be  combined  in  order  to  represent  a  complex  surface  efficiently  like 
the  one  of  a  ships  hull, 

For  the  particular  instance  of  the  harbor  tug  Huntington  (see  Fig.  13) 

the  more  advanced  methods,  which  have  been  applied,  are  a  modified 

version  of  the  "Gregory's  Square"  -  me t  h o d  that  can  even  interpolate 
drilateral  patches  as  a  superset  of  rectangular  patches  and  the 
■  method  for  arbitrary  triangles. 


CONCLUSI  ON 


The  current  version  of  the  two  phase  programming  system  has 
proven  that  the  program  organisation  can  meet  tne  expectations  for  a 
production  oriented  application.  The  objectives  of  portability, 
generality,  flexibility  and  efficiency  can  be  achieved. 

The  quality  of  the  interpolation  or,  better,  of  the 

representation  of  a  complex  shaped  surface  depends  very  much  on  the 
met  hods  chosen.  There  are  several  me t  h o d s  available  that  can  be 
incorporated  into  the  system  and  that  allow  to  combine  triangular  and 
rectangular  patches  and  thus  exploit  the  efficiency  of  the  rectangular 
interpolation  schemes  with  the  accuracy  of  global  shape  approximation 
ofthetriangularpatches. 

The  second  phase  of  the  current  version  of  the  programming  system 
was  designed  and  performed  for  the  purpose  of  illustration  rather  than 
for  final  application.  It  will  have  to  he  changed  for  a  particular 
application  of  graphical  or  production  device. 
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triangular  and  rectangular 
patches  (front  view) 


triangular  and  quadrilateral 
patches  (front  vi  ew) 


triangular  and  rectangular 
patches 

(front  view  from  the  inside) 


triangular  and  quadrilateral 
patches 

( f  ront  view  from  the  i  nsi  de) 


Fig.  13:  Results  from  improved  interpolation  methods  illustrated  on 
harbor  tug  Huntington 
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Since  the  output  data  of  the  first  phase  are  provided  in  the  most 
general  form  and  since  the  interpolation  points  are  computed  on  a 
regular  grid,  the  second  phase  of  the  programming  system  i  s  not 
limited  only  to  the  pure  interfacing  task.  It  also  could  include 
simple  numerical  integration  methods  in  order  to  provide  additional 
features,  such  as  the  computation  of  volumes,  centers  of  gravity  etc. 
of  various  parts  of  the  ship. 

The  system  needs  only  very  little  core  size  (10K  words) ,  so  that 
it  can  be  installed  even  on  smaller  computers  and  thus  offer  its 
benefits  also  to  the  small  shipbuilding  and  engineering  companies. 
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